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Abstract 

DAMA/LIBRA is running at the Gran Sasso National Laboratory of the 
I.N.F.N.. Here the results obtained with a further exposure of 0.34 ton x yr are 
presented. They refer to two further annual cycles collected one before and one af- 
ter the first DAMA/LIBRA upgrade occurred on September/October 2008. The 
cumulative exposure with those previously released by the former DAMA/Nal 
and by DAMA/LIBRA is now 1.17 ton x yr, corresponding to 13 annual cycles. 
The data further confirm the model independent evidence of the presence of 
Dark Matter (DM) particles in the galactic halo on the basis of the DM annual 
modulation signature (8.9 a C.L. for the cumulative exposure). In particular, 
with the cumulative exposure the modulation amplitude of the single-hit events 
in the (2 - 6) keV energy interval measured in Nal(Tl) target is (0.01 16 ±0.0013) 
cpd/kg/keV; the measured phase is (146 ± 7) days and the measured period is 
(0.999 ± 0.002) yr, values well in agreement with those expected for the DM 
particles. 



Keywords: Scintillation detectors, elementary particle processes, Dark Matter 

PACS numbers: 29.40.Mc - Scintillation detectors; 95.30.Cq - Elementary particle 
processes; 95.35.+d - Dark matter (stellar, interstellar, galactic, and cosmological) . 



1 Introduction 

The former DAMA/Nal [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14] and the present 
DAMA/LIBRA [15, 16, 17] experiments at the Gran Sasso National Laboratory have 
the main aim to investigate the presence of Dark Matter particles in the galactic 
halo by exploiting the model independent Dark Matter annual modulation signature 
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originally suggested in the mid 80's in ref. [18]. In fact, as a consequence of its annual 
revolution around the Sun, which is moving in the Galaxy travelling with respect to 
the Local Standard of Rest towards the star Vega near the constellation of Hercules, 
the Earth should be crossed by a larger flux of Dark Matter particles around ~ 2 
June (when the Earth orbital velocity is summed to the one of the solar system with 
respect to the Galaxy) and by a smaller one around ~ 2 December (when the two 
velocities are subtracted). Thus, this signature has a different origin and peculiarities 
than the seasons on the Earth and than effects correlated with seasons (consider the 
expected value of the phase as well as the other requirements listed below). This 
annual modulation signature is very distinctive since the effect induced by DM particles 
must simultaneously satisfy all the following requirements: the rate must contain a 
component modulated according to a cosine function (1) with one year period (2) and 
a phase that peaks roughly around ~ 2 nd June (3); this modulation must only be 
found in a well-defined low energy range, where DM particle induced events can be 
present (4); it must apply only to those events in which just one detector of many 
actually "fires" (single-hit events), since the DM particle multi- interaction probability 
is negligible (5); the modulation amplitude in the region of maximal sensitivity must 
be <7% for usually adopted halo distributions (6), but it can be larger in case of 
some possible scenarios such as e.g. those in refs. [19, 20]. This offers an efficient 
DM model independent signature, able to test a large interval of cross sections and of 
halo densities; moreover, the use of highly radiopure Nal(Tl) scintillators as target- 
detectors assures sensitivity to wide ranges of DM candidates, of interaction types and 
of astrophysical scenarios. 

It is worth noting that only systematic effects or side reactions able to simultane- 
ously fulfil all the 6 requirements given above (and no one has ever been suggested) 
and to account for the whole observed modulation amplitude might mimic this DM 
signature. 

The DAMA/LIBRA set-up, whose description, radiopurity and main features are 
discussed in details in ref. [16] has firstly been upgraded in September/October 2008: 
i) one detector has been recovered by replacing a broken PMT (see ref. [15]); ii) a new 
optimization of some PMTs and HVs has been performed; iii) all the transient digitizers 
recording the shape of the pulse have been replaced with new ones, the U1063A Acqiris 
8-bit lGS/s DC270 High-Speed cPCI Digitizers; iv) a new DAQ with optical read-out 
has been installed. Also during this upgrade the operations involving the handling of 
the sensitive part of the setup and the shield have been performed in HP Nitrogen 
atmosphere. The upgrade has allowed to enlarge the sensitive mass and to improve 
general features. Here we just remind that the sensitive part of this set-up is made of 
25 highly radiopure Nal(Tl) crystal scintillators (5-rows by 5-columns matrix) having 
9.70 kg mass each one. In each detector two 10 cm long special quartz light guides 
act also as optical windows on the two end faces of the crystal and are coupled to two 
low background photomultipliers working in coincidence at single photoelectron level. 
The detectors are housed in a sealed low-radioactive copper box installed in the center 
of a low-radioactive Cu/Pb/Cd-foils/polyethylene/paraffm shield; moreover, about 1 
m concrete (made from the Gran Sasso rock material) almost fully surrounds (mostly 
outside the barrack) this passive shield, acting as a further neutron moderator. A 
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threefold-levels sealing system excludes the detectors from the environmental air of 
the underground laboratory [16]. A hardware/software system to monitor the running 
conditions is operative and self-controlled computer processes automatically control 
several parameters and manage alarms. Moreover: i) the light response ranges from 
5.5 to 7.5 photoelectrons/keV, depending on the detector; ii) the hardware threshold 
of each PMT is at single photoelectron (each detector is equipped with two low back- 
ground photomultiplicrs working in coincidence); hi) energy calibration with X-rays/7 
sources are regularly carried out down to few keV; iv) the software energy threshold 
of the experiment is 2 keV; v) both single-hit events (where just one of the detectors 
fires) and multiple-hit events (where more than one detector fires) are acquired; v) the 
data are collected up to the MeV region despite the optimization is performed for the 
lower one. For the radiopurity, the procedures and further details see ref. [15, 16]. 

The data of the former DAMA/Nal (0.29 ton x yr) and those of the first 4 annual 
cycles of DAMA/LIBRA (total exposure 0.53 tonxyr) have already given positive 
model independent evidence for the presence of DM particles in the galactic halo with 
high confidence level on the basis of the DM annual modulation signature [15]. 

In this paper the model independent results with other two annual cycles 
DAMA/LIBRA-5,6 are presented. As mentioned, the data of the first cycle have 
been collected in the same conditions as DAMA/LIBRA-1,2,3,4 [15, 16], while the 
data of DAMA/LIBRA-6 have been taken after the above mentioned 2008 upgrade. 

2 The results 

The updated exposures of the DAMA/LIBRA annual cycles and the cumulative one 
with the former DAMA/Nal are given in Table 1. 

Table 1: Exposures of the DAMA/LIBRA-5,6 annual cycles. Here a = (cos 2 w(t — t )) 
is the mean value of the squared cosine and (3 — (cosuj(t — 1 )} is the mean value of the 
cosine (the averages are taken over the live time of the data taking and to — 152.5 day, 
i.e. June 2 nd ); thus, (a — (3 2 ) indicates the variance of the cosine (i.e. it is 0.5 for a full 
year of data taking). The information on the previously published DAMA/LIBRA- 
1,2,3,4, are recalled as well as the cumulative exposure, when including the former 
DAMA/Nal. 





Period 


mass 
(kg) 


Exposure 
(kg x day) 




DAMA/LIBRA-1,2,3,4 


Sept. 9, 2003 - July 17, 2007 


232.8 


192824 


0.537 


DAMA/LIBRA-5 


July 17, 2007 - Aug. 29, 2008 


232.8 


66105 


0.468 


DAMA/LIBRA-6 


Nov. 12,2008 - Sept, 1, 2009 


242.5 


58768 


0.519 


DAMA/LIBRA- 1 to -6 


Sept. 9, 2003 - Sept. 1, 2009 




317697- 0.87 tonxyr 


0.519 


DAMA/Nal + DAMA/LIBRA- 1 to 6: 




1.17 tonxyr 





The only data treatment, which is performed on the raw data, is to remove noise 
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pulses (mainly PMT noise, Cherenkov light in the light guides and in the PMT win- 
dows, and afterglows) near the energy threshold in the single-hit events; for a descrip- 
tion of the used procedure and details sec ref. [16]. 

In the DAMA/LIBRA-1,2,3,4,5,6 annual cycles about 7.2 x 10 7 events have also 
been collected for energy calibrations and about 3 x 10 6 events/keV for the evaluation 
of the acceptance windows efficiency for noise rejection near energy threshold. The 
periodical calibrations and, in particular, those related with the acceptance windows 
efficiency mainly affect the duty cycle of the experiment. From Table 1 one can observe 
a significant improvement in the duty cycle of the sixth annual cycles with respect to 
the previous ones; this is mainly due to the new transient Digitizers and DAQ installed 
at fall 2008. 

Several analyses on the model-independent investigation of the DM annual modu- 
lation signature have been performed as previously done in ref. [15] and refs. therein. 
In particular, Fig. 1 shows the time behaviour of the experimental residual rates for 
single-hit events in the (2-4), (2-5) and (2-6) keV energy intervals. These residual 
rates are calculated from the measured rate of the single-hit events (already corrected 
for the overall efficiency and for the acquisition dead time) after subtracting the con- 
stant part: < r^-fe — flatjk >jk- Here is the rate in the considered i-th time 
interval for the j-th detector in the fc-th energy bin, while flatjk is the rate of the 
j-th detector in the k-th energy bin averaged over the cycles. The average is made 
on all the detectors (j index) and on all the energy bins (k index) which constitute 
the considered energy interval. The weighted mean of the residuals must obviously be 
zero over one cycle. 

For clarity in Fig. 1 only the DAM A/LIBRA data collected over six annual cycles 
(0.87 ton x yr) are shown; the DAMA/Nal data (0.29 ton x yr) and comparison with 
DAMA/LIBRA are available in ref. [15]. 

The hypothesis of absence of modulation in the data can be discarded (see Table 

2). 

Table 2: Test of absence of modulation in the data of the DAMA/LIBRA-1,2,3,4,5,6 
and without /with also the data of the former DAMA/Nal. As it can be seen, a null 
modulation amplitude is discarded by the data. 



Energy interval 
(keV) 


DAMA/LIBRA 
(6 annual cycles) 


DAMA/Nal & DAMA/LIBRA 
(7+6 annual cycles) 


2-4 


X 2 /d.o.f. = 90.0/43 
->■ P = 3.6 x 10" 5 


xVd.o.f. = 147.4/80 
-> P = 6.8 x 10~ 6 


2-5 


X 2 /d.o.f. = 82.1/43 
->■ P = 3.1 xlO" 4 


X 2 /d.o.f. = 135.2/80 
— >• P = 1.1 x 10~ 4 


2-6 


X 2 /d.o.f. = 68.9/43 
->■ P = 7.4 x 10~ 3 


xVd.o.f. = 139.5/80 
-> P = 4.3 x 10~ 5 



The single-hit residual rate of DAMA/LIBRA-1,2,3,4,5,6 of Fig. 1 can be fitted 
with the formula: Acosoj(t — to) considering a period T = — = 1 yr and a phase 
to = 152.5 day (June 2 nd ), as expected by the DM annual modulation signature; this 
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Figure 1: Experimental model-independent residual rate of the single-hit scintillation 
events, measured by DAMA/LIBRA,1,2,3,4,5,6 in the (2 - 4), (2 - 5) and (2 - 6) 
keV energy intervals as a function of the time. The zero of the time scale is January 
I s * of the first year of data taking of the former DAMA/Nal experiment [15]. The 
experimental points present the errors as vertical bars and the associated time bin 
width as horizontal bars. The superimposed curves are the cosinusoidal functions 
behaviors Acosoo(t — to) with a period T = — = 1 yr, with a phase to — 152.5 day 
(June 2 nd ) and with modulation amplitudes, A, equal to the central values obtained 
by best fit over the whole data including also the exposure previously collected by 
the former DAMA/Nal experiment: cumulative exposure is 1.17 ton x yr (see also 
ref. [15] and refs. therein). The dashed vertical lines correspond to the maximum 
expected for the DM signal (June 2 nd ), while the dotted vertical lines correspond to 
the minimum. See text. 
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can be repeated for the total available exposure 1.17 ton x yr including the former 
DAMA/Nal data (see [15] and refs. therein). The results are shown in Table 3. 

Tabic 3: Modulation amplitude, A, obtained by fitting the single-hit residual rate of the 
six DAM A/LIBRA annual cycles (Fig. 1), and including also the former DAMA/Nal 
data given elsewhere (see [15] and refs. therein) for a total cumulative exposure of 1.17 
ton x yr. It has been obtained by fitting the data with the formula: Acosco(t — to) 
with T = 22: = 1 yr and to — 152.5 day (June 2 nd ), as expected for a signal by the 
DM annual modulation signature. The corresponding x 2 value for each fit and the 
confidence level are also reported 



Energy interval 
(keV) 


DAM A/LIBRA 

(cpd/kg/keV) 


DAMA/Nal & DAM A/LIBRA 

(cpd/kg/keV) 


2-4 


A=(0.0170±0.0024) 
X 2 /d.o.f. = 41.0/42 


A=(0.0183±0.0022) 
X 2 /d-o.f. = 75.7/79 
-> 8.3 a C.L. 


2-5 


A=(0.0129±0.0018) 
X 2 /d-o.f. = 30.7/42 


A=(0.0144±0.0016) 
X 2 /d.o.f. = 56.6/79 
— > 9.0 cr C.L. 


2-6 


A=(0.0097±0.0015) 
X 2 /d-o.f. = 24.1/42 


A=(0.0114±0.0013) 
X 2 /d.o.f. = 64.7/79 
— > 8.8 cr C.L. 



The compatibility among the 13 annual cycles has been investigated. In particular, 
the modulation amplitudes measured in each annual cycle of the whole 1.17 ton x yr 
exposure have been analysed as in ref. [15]. Indeed these modulation amplitudes are 
normally distributed around their best fit value as pointed out by the x 2 test (x 2 = 9.3, 
12.2 and 10.1 over 12 d.o.f. for the three energy intervals, respectively) and the run 
test (lower tail probabilities of 57%, 47% and 35% for the three energy intervals, re- 
spectively). Moreover, the DAMA/LIBRA-5 and DAM A/LIB RA-6 (2-6) keV modu- 
lation amplitudes are (0.0086±0.0032) cpd/kg/keV and (0.0101±0.0031) cpd/kg/kcV, 
respectively, in agreement with that of DAMA/LIBRA-1,2,3,4: (0.0110 ± 0.0019) 
cpd/kg/keV; we also recall that the statistical compatibility between the DAMA/Nal 
and DAMA/LIBRA-1,2,3,4 modulation amplitudes has been verified [15]. Thus, 
also when adding DAMA/LIBRA-5, 6, the cumulative result from DAMA/Nal and 
DAMA/LIBRA can be adopted. 

Table 4 shows the results obtained for the cumulative 1.17 ton x yr exposure when 
the period and phase parameters are kept free in the fitting procedure described above. 

The DAMA/LIBRA single-hit residuals of Fig.l and those of DAMA/Nal (see 
e.g. [15]) have also been investigated by a Fourier analysis, obtaining a clear peak 
corresponding to a period of 1 year (see Fig. 2); the same analysis in other energy 
region shows instead only aliasing peaks. 

The measured energy distribution has been investigated in other energy regions 
not of interest for Dark Matter, also verifying the absence of any significant back- 
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Table 4: Modulation amplitude (A), period (T = — ) and phase (to), obtained by 
fitting, with the formula: A cos ui(t — to), the single-hit residual rate of the cumulative 
1.17 ton x yr exposure. The results are well compatible with expectations for a signal 
in the DM annual modulation signature. 



Energy interval 


A (cpd/kg/keV) 


T = 77 (yr) 


t (days) 


C. L. 


2-4 


(0.0194±0.0022) 


(0.996±0.002) 


136±7 


8.8ct 


2-5 


(0.0149±0.0016) 


(0.997±0.002) 


142±7 


9.3cr 


2-6 


(0.0116±0.0013) 


(0.999±0.002) 


146±7 


8.9ct 




Frequency (d" 1 ) Frequency (d" 1 ) 



Figure 2: Power spectrum of the measured single-hit residuals in the (2-6) keV (solid 
lines) and (6-14) keV (dotted lines) energy intervals calculated according to ref. [21], 
including also the treatment of the experimental errors and of the time binning. The 
data refer to: a) DAMA/LIBRA-1,2,3,4,5,6 (exposure of 0.87 ton x yr); b) the cumu- 
lative 1.17 ton x yr exposure (DAMA/Nal and DAMA/LIBRA-1,2,3,4,5,6). As it can 
be seen, the principal mode present in the (2-6) keV energy interval corresponds to 
a frequency of 2.697 x 10 -3 d~ 1 and 2.735 x 10~ 3 d _1 (vertical lines), respectively in 
the a) and b) case. They correspond to a period of ~ 1 year. A similar peak is not 
present in the (6-14) keV energy interval just above. 

ground modulation 1 . Following the procedures described in ref. [15] and ref. therein, 
the measured rate integrated above 90 keV, R90, as a function of the time has been 
analysed. In particular, also for these two latter annual cycles the distribution of the 
percentage variations of R90 with respect to the mean values for all the detectors has 
been considered; it shows a cumulative gaussian behaviour with er ~ 1%, well ac- 



1 , In fact, the background in the lowest energy region is essentially due to "Compton" electrons, 
X-rays and/or Auger electrons, muon induced events, etc., which are strictly correlated with the 
events in the higher energy part of the spectrum. Thus, if a modulation detected in the lowest energy 
region would be due to a modulation of the background (rather than to a signal) , an equal or larger 
modulation in the higher energy regions should be present. 
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counted by the statistical spread expected from the used sampling time (see Fig. 3). 
Moreover, fitting the time behaviour of R90 with phase and period as for DM particles, 
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Figure 3: Distribution of the percentage variations of R90 with respect to the mean 
values for all the detectors in the DAMA/LIBRA-5,6 annual cycles (histogram); the 
superimposed curve is a gaussian fit. 

a modulation amplitude compatible with zero is also found in DAMA/LIBRA-5 and 
DAMA/LIBRA-6: (0.20 ±0.18) cpd/kg and (-0.20 ±0.16) cpd/kg, respectively. This 
also excludes the presence of any background modulation in the whole energy spec- 
trum at a level much lower than the effect found in the lowest energy region for the 
single-hit events. In fact, otherwise - considering the R90 mean values - a modulation 
amplitude of order of tens cpd/kg, that is ~ 100 a far away from the measured value, 
would be present. Similar result is obtained when comparing the single-hit residuals 
in the (2-6) keV with those in other energy intervals; see as an example Fig. 4. It 
is worth noting that the obtained results already account for whatever kind of back- 
ground and, in addition, that no background process able to mimic the DM annual 
modulation signature (that is able to simultaneously satisfy all the peculiarities of the 
signature and to account for the measured modulation amplitude) is available (see also 
discussions e.g. in [15, 22]). 

A further relevant investigation has been performed by applying the same hardware 
and software procedures, used to acquire and to analyse the single-hit residual rate, 
to the multiple-hit one. In fact, since the probability that a DM particle interacts in 
more than one detector is negligible, a DM signal can be present just in the single-hit 
residual rate. Thus, the comparison of the results of the single-hit events with those 
of the multiple-hit ones corresponds practically to compare between them the cases of 
DM particles beam-on and beam-off. This procedure also allows an additional test 
of the background behaviour in the same energy interval where the positive effect is 
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Figure 4: Experimental residuals in the (2 - 6) keV region and those in the (6 - 14) 
keV energy region just above for the cumulative 1.17 ton x yr, considered as collected 
in a single annual cycle. The experimental points present the errors as vertical bars 
and the associated time bin width as horizontal bars. The initial time of the figure is 
taken at August 7 th . The clear modulation satisfying all the peculiarities of the DM 
annual modulation signature is present in the lowest energy interval, while it is absent 
just above; in fact, in the latter case the best fitted modulation amplitude is: (0.00007 
± 0.00077) cpd/kg/keV. 



observed. In particular, in Fig. 5 the residual rates of the single-hit events measured 
over the six DAMA/LIBRA annual cycles are reported, as collected in a single cycle, 
together with the residual rates of the multiple-hit events, in the considered energy 
intervals. While, as already observed, a clear modulation, satisfying all the peculiarities 
of the DM annual modulation signature, is present in the single-hit events, the fitted 
modulation amplitudes for the multiple-hit residual rate are well compatible with zero: 
(-0.0011 ± 0.0007) cpd/kg/keV, (-0.0008 ± 0.0005) cpd/kg/keV, and (-0.0006 ± 
0.0004) cpd/kg/keV in the energy regions (2 - 4), (2 - 5) and (2 - 6) keV, respectively. 
Thus, again evidence of annual modulation with proper features as required by the 
DM annual modulation signature is present in the single-hit residuals (events class to 
which the DM particle induced events belong), while it is absent in the multiple-hit 
residual rate (event class to which only background events belong). Similar results 
were also obtained for the last two annual cycles of the DAMA/Nal experiment [6]. 
Since the same identical hardware and the same identical software procedures have 
been used to analyse the two classes of events, the obtained result offers an additional 
strong support for the presence of a DM particle component in the galactic halo. 

As in ref. [15], the annual modulation present at low energy can also be shown by 
depicting - as a function of the energy - the modulation amplitude, S m ,k, obtained 
by maximum likelihood method over the data considering T —1 yr and to = 152.5 
day. For such purpose the likelihood function of the single-hit experimental data in 

the k— th energy bin is defined as: Lk = Ti^e MiJ ' fc I ^' k k , , where Nijk is the number of 
events collected in the i-th time interval (hereafter 1 day), by the j-th detector and 
in the k-th energy bin. N^k follows a Poisson's distribution with expectation value 
fJ-ijk = [bjk + Sik] MjAtiAEejk- The bjk are the background contributions, Mj is 
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Figure 5: Experimental residual rates over the six DAM A/LIBRA annual cycles for 
single-hit events (open circles) (class of events to which DM events belong) and for 
multiple-hit events (filled triangles) (class of events to which DM events do not belong). 
They have been obtained by considering for each class of events the data as collected 
in a single annual cycle and by using in both cases the same identical hardware and the 
same identical software procedures. The initial time of the figure is taken on August 
7 th . The experimental points present the errors as vertical bars and the associated 
time bin width as horizontal bars. See text and ref. [15]. Analogous results were 
obtained for the DAMA/Nal data [6]. 



the mass of the j— th detector, At, is the detector running time during the i-th time 
interval, AE is the chosen energy bin, €jk is the overall efficiency. Moreover, the signal 
can be written as Sik = So,k + Sm.fc ■ cosu(ti — to), where So,fc is the constant part of 
the signal and S m .k is the modulation amplitude. The usual procedure is to minimize 
the function = — 2^n(Lk) — const for each energy bin; the free parameters of the fit 
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are the (bjk + So,k) contributions and the S m ,k parameter. Hereafter, the index k is 
omitted when unnecessary. 
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Figure 6: Energy distribution of the S m variable for the total cumulative exposure 
1.17 tonxyr. The energy bin is 0.5 keV. A clear modulation is present in the lowest 
energy region, while S m values compatible with zero are present just above. In fact, 
the S m values in the (6-20) keV energy interval have random fluctuations around zero 
with x 2 equal to 27.5 for 28 degrees of freedom. 

In Fig. 6 the obtained S m are shown in each considered energy bin (there AE = 0.5 
keV). It can be inferred that positive signal is present in the (2-6) keV energy interval, 
while S m values compatible with zero are present just above. In fact, the S m values 
in the (6-20) keV energy interval have random fluctuations around zero with x 2 equal 
to 27.5 for 28 degrees of freedom. All this confirms the previous analyses. 



The method also allows the extraction of the the S m values for each detector, for 
each annual cycle and for each energy bin. Thus, following the procedure described 
in ref. [15], we have also verified that the S TO are statistically well distributed in all 
the six DAMA/LIBRA annual cycles and in all the sixteen energy bins (AE = 0.25 
keV in the 2-6 keV energy interval) for each detector. Moreover, that procedure also 
allows the definition of a x 2 f° r each detector; the associated degree of freedom are 
16 for the detector restored after the upgrade in 2008 and 96 for the others. The 
values of the x 2 /d.o.f. range between 0.7 and 1.22 for twenty-four detectors, and the 
observed annual modulation effect is well distributed in all these detectors at 95% 
C.L.. A particular mention is deserved to the remaining detector whose value is 1.28 
exceeding the value corresponding to that C.L.; this also is statistically consistent, 
considering that the expected number of detector exceeding this value over twenty- 
five is 1.25. Moreover, the mean value of the 25 x 2 /d-o.f. is 1.066, slightly larger 
than expected. Although this can be still ascribed to statistical fluctuations (see 
before), let us ascribe it to a possible systematics. In this case, one would have an 
additional error of < 4 x 10~ 4 cpd/kg/keV, if quadratically combined, or < 5 x 10~ 5 
cpd/kg/keV, if linearly combined, to the modulation amplitude measured in the (2 - 
6) keV energy interval. This possible additional error: < 4% or < 0.5%, respectively, 
of the DAMA/LIBRA modulation amplitude is an upper limit of possible systematic 
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effects. 

Among further additional tests, the analysis of the modulation amplitudes as a 
function of the energy separately for the nine inner detectors and the remaining exter- 
nal ones has been carried out including the DAMA/LIBRA-5,6 data to those already 
analysed in ref. [15]. The obtained values are fully in agreement; in fact, the hypoth- 
esis that the two sets of modulation amplitudes as a function of the energy belong to 
same distribution has been verified by x 2 test, obtaining: x 2 /d.o.f. — 3.1/4 and 7.1/8 
for the energy intervals (2-4) and (2-6) keV, respectively (AE = 0.5 keV). This shows 
that the effect is also well shared between inner and external detectors. 

Let us, finally, release the assumption of a phase to = 152.5 day in the procedure 
to evaluate the modulation amplitudes from the data of the 1.17 ton x yr. In this case 
alternatively the signal has been written as: 

Sik = #o,fc + £m,fcCosw(fj-t Q ) + Z TO fcsmo^tj-io) = S k + Y m k cosu}(ti-t*),(l) 



For signals induced by DM particles one would expect: i) Z m> u ~ (because of 
the orthogonality between the cosine and the sine functions); ii) S m ,k — Y m ,k> iii) 
t* ~ to = 152.5 day. In fact, these conditions hold for most of the dark halo models; 
however, it is worth noting that slight differences can be expected in case of possible 
contributions from non-thermalized DM components, such as e.g. the SagDEG stream 
[8] and the caustics [23]. 
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Figure 7: 2a contours in the plane (S m ,Z m ) (left) and in the plane (Y m ,t*) (right) 
for the (2-6) keV and (6-14) keV energy intervals. The contours have been obtained 
by the maximum likelihood method, considering the cumulative exposure of 1.17 ton 
x yr. A modulation amplitude is present in the lower energy intervals and the phase 
agrees with that expected for DM induced signals. 



Fig. 7 -left shows the 2a contours in the plane (S mi Z m ) for the (2-6) keV and (6- 
14) keV energy intervals and Fig. 7 -right shows, instead, those in the plane (Y m ,t*). 
Table 5 shows the best fit values for the (2-6) and (6-14) keV energy interval (la 
errors) for S m versus Z m and Y m versus t*. 

Finally, forcing to zero the contribution of the cosine function in eq. (1), the Z m 
values as function of the energy have also been determined by using the same procedure. 
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Tabic 5: Best fit values for the (2-6) and (6-14) keV energy interval (la errors) for 
S m versus Z m and Y m versus t*, considering the cumulative exposure of 1.17 ton x 
yr. See also Fig. 7. 



E 
(keV) 


Sm 

(cpd/kg/keV) 


Z m 

(cpd/kg/keV) 


Y m 

(cpd/kg/keV) 


t* 
(day) 


2-6 
6-14 


(0.0111 ± 0.0013) 
-(0.0001 ± 0.0008) 


-(0.0004 ± 0.0014) 
(0.0002 ± 0.0005) 


(0.0111 ± 0.0013) 
-(0.0001 ± 0.0008) 


(150.5 ± 7.0) 
undefined 



The values of Z m as a function of the energy is reported in Fig. 8. Obviously, such 
values are expected to be zero in case of presence of a DM signal with t* ~ to = 152.5 
day. By the fact, the % 2 test applied to the data supports the hypothesis that the Z m 
values are simply fluctuating around zero; in fact, for example in the (2-14) keV and 
(2-20) keV energy region the x 2 /d.o.f. are equal to 21.6/24 and 47.1/36 (probability 
of 60% and 10%), respectively. 
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Figure 8: Energy distribution of the Z m variable for the total exposure (1.17 ton x 
yr, DAM A/Nal&DAM A/LIBRA), once forced to zero the contribution of the cosine 
function in eq. (1). The energy bin is 0.5 keV. The Z m values are expected to be zero 
in case of presence of a DM particles' signal with t* ~ to — 152.5 day. By the fact, 
the x 2 test applied to the data supports the hypothesis that the Z m values are simply 
fluctuating around zero; see text. 



The behaviours of the Y m and of the phase t* variables as function of energy are 
shown in Fig. 9 for the total exposure (1.17 ton x yr, D AM A/Nal&D AM A/LIBRA). 
The Y m arc superimposed with the S m values with 1 keV energy bin (unlike Fig. 6 
where the energy bin is 0.5 keV). As in the previous analyses, an annual modulation 
effect is present in the lower energy intervals and the phase agrees with that expected 
for DM induced signals. 

These results confirm those achieved by other kinds of analyses. 

Sometimes naive statements were put forwards as the fact that in nature several 
phenomena may show some kind of periodicity. It is worth noting that the point is 
whether they might mimic the annual modulation signature in DAM A/LIBRA (and 
former DAMA/Nal), i.e. whether they might be not only quantitatively able to account 
for the observed modulation amplitude but also able to contemporaneously satisfy all 
the requirements of the DM annual modulation signature. The same is also for side 
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Figure 9: Top: Energy distributions of the Y m variable (light data points; red colour 
online) and of the S m variable (solid data points) for the total exposure (1.17 ton x 
yr, DAMA/Nal&DAMA/LIBRA). Here, unlike the data of Fig. 6, the energy bin is 
1 keV. Bottom: Energy distribution of the phase t* for the total exposure; here the 
errors are at 2a. An annual modulation effect is present in the lower energy intervals 
up to 6 keV and the phase agrees with that expected for DM induced signals. No 
modulation is present above 6 keV and the phase is undetermined. 



reactions. This has already been deeply investigated in ref. [15, 16] and references 
therein; the arguments and the quantitative conclusions, presented there, also apply 
to the DAMA/LIBRA-5,6 data. Some additional arguments have also been recently 
addressed in [22, 24]. 



3 Comments 

The obtained model independent evidence - at 8.9 a C.L. over 13 annual cycles - is 
compatible with a wide set of scenarios regarding the nature of the DM candidate and 
related astrophysical, nuclear and particle Physics (see e.g. ref. [5, 6, 7, 9, 10, 11, 12], 
Appendix A of ref. [15] and in literature, for example see [19, 25, 26]); and many other 
possibilities are open. Further future works are foreseen. 

It is worth recalling that no other experiment exists, whose result can be directly 
compared in a mo del- independent way with those by DAMA/Nal and DAMA/LIBRA, 
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and that - more in general - results obtained with different target materials and/or 
different approaches cannot be directly compared among them in a model-independent 
way. This is in particular due to the existing experimental and theoretical uncertain- 
ties, not last e.g. how many kinds of dark matter particles can exist in the Universe 2 , 
the nature, the interaction types, the different nuclear and/or atomic correlated as- 
pects, the unknown right halo model, the right DM density, etc. as well as the uncer- 
tainties on the values of each one of the many involved experimental and theoretical 
parameter/assumption/approximation used in the calculations. Moreover, some ex- 
perimental aspects of some techniques used in the field have also to be addressed 
[5, 27, 24]. Another relevant argument is the methodological robustness [28]. In par- 
ticular, the general considerations on comparisons reported in Appendix A of ref. [15] 
still hold. Hence, claims for contradiction have no scientific basis. On the other hand, 
whatever possible "positive" result has to be interpreted and a large room of compat- 
ibility with DAMA annual modulation evidence is present. 

Similar considerations can also be done for the indirect detection searches, since 
it does not exist a biunivocal correspondence between the observables in the direct 
and indirect experiments. However, if possible excesses in the positron to electron flux 
ratio and in the 7 rays flux with respect to a modeling of the background contribution, 
which is expected from the considered sources, might be interpreted - under some 
assumptions - in terms of Dark Matter, this would also be not in conflict with the effect 
observed by DAMA experiments. It is worth noting that different possibilities either 
considering different background modeling or accounting for other kinds of sources can 
also explain the indirect observations [29] . 

Finally, as regards the accelerator searches for new particles beyond the Standard 
Model of particle Physics, it is worth noting that they can demonstrate the existence 
of some of the possible DM candidates, but cannot credit that a certain particle is the 
DM solution or the "single" DM solution. Moreover, DM candidates and scenarios 
exist (even e.g. for the neutralino candidate) on which accelerators cannot give any 
information. It is also worth noting that for every candidate (including the neutralino) 
there exist various different possibilities for the theoretical aspects. Nevertheless, the 
results from accelerators will give outstanding and crucial complementary information 
in the field. 

A new upgrade of DAMA/LIBRA is foreseen in 2010 with the replacement of all 
the low background PMTs with new ones having higher quantum efficiency; the main 
aim is to lower the software energy threshold and, thus, to increase the experimental 
sensitivity and to disentangle - in the corollary investigation on the candidate parti- 
cle^) - at least some of the many possible astrophysical, nuclear and particle Physics 
scenarios and related experimental and theoretical uncertainties. 



2 In fact, it is worth noting that, considering the richness in particles of the visible matter which is 
less than 1% of the Universe density, one could also expect that the particle part of the Dark Matter 
in the Universe may also be multicomponcnt. 
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4 Conclusions 



The new annual cycles DAMA/LIBRA-5,6 have further confirmed a peculiar annual 
modulation of the single-hit events in the (2-6) keV energy region satisfying the 
many requests of the DM annual modulation signature; the total exposure by for- 
mer DAMA/Nal and present DAMA/LIBRA is 1.17 ton x yr. 

In fact, as required by the DM annual modulation signature: 1) the single-hit events 
show a clear cosine-like modulation as expected for the DM signal; 2) the measured 
period is equal to (0.999 ± 0.002) yr well compatible with the 1 yr period as expected 
for the DM signal; 3) the measured phase (146 ± 7) days is well compatible with the 
roughly ~ 152.5 days expected for the DM signal; 4) the modulation is present only 
in the low energy (2-6) keV energy interval and not in other higher energy regions, 
consistently with expectation for the DM signal; 5) the modulation is present only in 
the single-hit events, while it is absent in the multiple-hit ones as expected for the DM 
signal; 6) the measured modulation amplitude in Nal(Tl) of the single-hit events in the 
(2-6) keV energy interval is: (0.0116 ± 0.0013) cpd/kg/keV (8.9 a C.L.). No system- 
atic or side processes able to simultaneously satisfy all the many peculiarities of the 
signature and to account for the whole measured modulation amplitude is available. 
Further work is in progress. 
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